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ABSTRACT 

We present and analyze high-resolution optical spectra of two Algol binaries, 
namely S Equ and KO Aql, obtained with the echelle spectrograph at Catania As- 
trophysical Observatory. New accurate radial velocities for the hotter primary com- 
l/-) ' ponents are obtained. Thanks to the cross-correlation procedure, we were able to 
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measure, for the first time to our knowledge, radial velocities also for the cool sec- 



ondary components of S Equ and KO Aql. By combining the parameters obtained 
from the solution of the radial velocity curves with those obtained from the light 
J""""- | curve analysis, reliable absolute parameters of the systems have been derived. The 



rotational velocity of the hotter components of S Equ and KO Aql has been measured 
and it is found that the gainers of both systems rotate about 30% faster than syn- 
chronously. This is likely due to mass transfer across the Lagrangian L\ point from the 
cooler to the hotter component. The lower luminosity of the mass- gaining components 
of these systems compared to normal main-sequence stars of the same mass can be 
also an effect of the mass transfer. The Ha profiles were analyzed with the "synthesis 
and subtraction" technique and reveal clear evidence of mass transfer and accretion 
structures. In both systems, especially before the primary eclipses and afterwards, we 
clearly observed extra-absorption lines. From the integrated absorption and the radial 
velocity variations of these features, we found that the mass accretion is very dense 
around the impact region of the hotter components. A double-peaked emission in the 
spectra of S Equ was seen outside the eclipses. One of these peaks is likely originated 
in a region between the center of mass and the cooler component, which is occupied 
by the flowing matter. Furthermore, the Ha difference spectra of S Equ and KO Aql 
display also emission features, which should be arising from the magnetic activity of 
the cooler components. 

Key words: stars: binaries: eclipsing - stars: mass transfer - stars: individual: S Equ, 
KO Aql 
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1 INTRODUCTION 

Algol systems are close, interacting binaries, which consist 
of a cool, giant or subgiant secondary star of F-K spectral 
type that fills in its own Roche lobe and transfers mass and 
angular momentum to a hot and more massive B-A type 
main-sequence primary star. They are an important astro- 
physical laboratory for the study of different phenomena 
such as mass transfer and accretion, magnetic activity in 
the late-type companion, angular momentum and orbit evo- 
lution. A crucial need for the study of all these processes is 
the knowledge of the basic stellar and orbital parameters. In 
the last decades, the number of Algol-type binaries with reli- 
able fundamental parameters has considerably increased due 
to the availability of high-resolution spectroscopic data an- 
alyzed with new powerful techniques (e.g. cross-correlation, 
Popper & Jeong, 1994) to measure radial velocities, and to 
the accurate analysis of multi-band (from visible to near- 
IR wavelength) light curves with sophisticated spot models 
(e.g., Lazaro et al. 2004). 

Observational evidences of mass transfer activity and 
accretion structures in Algols are: (i) the variations of the 
light and color curves, (ii) distortions of the radial velocity 
curves of the primary components, (iii) the mass-gainer hot- 
ter components rotating faster than synchronously, and (iv) 
extra emission and/or absorption features in strong spectral 
lines such as Ha and H/3. Some theoretical works have tack- 
led the problem of Roche overflow and mass transfer, pre- 
dicting the path of the gas stream and the accretion struc- 
tures around the primary component (see, e.g., Prendergast 
& Taam 1974, Lubow & Shu 1975, Peters 1989). 

In more recent years, different spectral behaviors (espe- 
cially extra Ha emission and/or absorption features) have 
been observed for short- and long-period Algol systems and 
have been attributed to different effects of mass transfer 
due to its strength and to the system geometry (Richards 
& Albright 1999). In long-period systems (P > 6 days), 
the gas stream misses the primary component and give rise 
to a "permanent accretion disk", because the size of the 
mass gainer is very small relative to the binary separation. 
The mass transfer and accretion process is more complex in 
short-period Algols (P ^ 6 days). In this case, the stream 



particles mostly impact the primary's photosphere due to 
the large size of the primary relative to the separation of the 
components. The accretion structure is called a "transient 
accretion disk" in Algols with orbital periods in the range 
4.5-6 days (Kaitchuck & Honeycutt 1982) and an "accretion 
annulus" in the systems with P or b < 4.5 days (Richards et 
al., 1996). A wide spectroscopic survey of short- and long- 
period Algol systems aimed at the study of mass transfer has 
been performed by Richards & Albright (1999) and Vesper 
et al. (2001). The two Algol systems studied by us in the 
present work are both short-period ones. 

S Equulei (BD +04°4584, HD 199454, HIP 103419, V~ 
8™4) is a semidetached Algol-type eclipsing binary with an 
orbital period of P = 3.4361 days. The hotter and cooler 
components have been classified as B9.5 V and F9 III-IV, 
respectively, by Plavec (1966) and Plavec & Polidan (1976). 

The first complete light curve was presented by Cata- 
lano & Rodono (1968) who also obtained the light curve so- 
lution. They found the system to be partially eclipsing and 
observed some peculiarities in the light curve. In particular, 
they noticed a decrease in luminosity, especially in B light 
(about 0™04), occurring near 0.9 phase. They interpreted 
this anomaly as due to the gaseous stream. Piotrowski et 
al. (1974) have discussed the distortions of the light curves 
of some semi-detached systems, including S Equ with the 
Catalano & Rodono data. They ascribed the brightening 
observed in several systems just after the primary minimum 
to emission from the gas stream that is observed along its 
axis, in the direction of higher optical depth. The luminosity 
decrease observed before the primary minimum is attributed 
to absorption by electron scattering produced by the stream 
itself that is projected against the disk of the hotter com- 
ponent at these phases. A more recent light curve solution 
(Zola, 1992) suggests a third light comparable to the con- 
tribution of the cooler component to the total flux of the 
system. 

Richards & Albright (1999), Vesper et al. (2001), and 
Richards (2001) have studied S Equ spectroscopically and 
found single-peaked Ha emission features with additional 
weak double-peaked emission at some epochs, indicating 
mass transfer activity and a transient disk structure around 
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the primary component. Doppler tomography based on 
the Ha differences spectra displays extra-absorption and/or 
emission due to the gas stream, accretion structures, and 
chromospheric activity associated with the cooler compo- 
nent. 

A detailed O-C analysis of the times of minima of S Equ 
was recently performed by Soydugan et al. (2003). They 
found an increase in the orbital period of the system with a 
rate of 0.102 syr -1 , requiring a mass transfer rate dM/dt ~ 
3.97 x 10 -8 Mq yr" 1 from the secondary to the primary 
component. 

KO Aquilae (BD +10°3655, HD 92177, V~ 8 m 4) is also 
a classical Algol-type binary. Its variability was discovered 
by Hoffmeister (1930) and its orbital period was determined 
as 2.864071 days by Plaut (1932). Sahade (1945) estimated 
the spectral type of the primary component as AO. Blanco 
& Cristaldi (1974) and Mader & Angione (1996) analyzed 
the light curves of the system and determined the stellar 
parameters. In the latter study, an orbital period increase 
with a rate of 0.375 syr -1 was also determined and explained 
as due to the mass transfer process. Sahade (1945) published 
the radial velocity curve and orbit solution of the primary 
component using plate spectra. He was not able to detect 
any sign of the cool component in his spectra. Vesper et 
al. (2001) obtained only one spectrum in the Ha region, 
at the orbital phase 0.611, and observed extra absorption 
and emission attributed to accretion structures around the 
primary component. 

The basic orbital and physical parameters of S Equ and 
KO Aql, as compiled from the literature, are reported in 
Table 1. 

In this work, we present the results of the analysis of 
new high-resolution echelle spectra of S Equ and KO Aql. 
The observations and data reduction are described in §2. 
We report the radial velocity measurements and the first 
complete orbit solutions for the components of S Equ and 
KO Aql in §3. The projected rotational velocities of the hot- 
ter components are discussed in §4. The results of the ap- 
plication of synthetic spectra to reproduce the photospheric 
contribution of both components, to refine the spectral clas- 
sification, and to study the Ha extra-absorption and emis- 
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Table 1. Basic parameters of S Equ and KO Aql 



System 


Spectral 


P orb 


i 


Ti 


T 2 


ri 


T2 


Ref. 




Type 


(days) 


(deg) 


(K) 


(K) 


(=Ri/a) 


(=R 2 /a) 




S Equ 


B8-9.5V -j 


- F9III-IV 3.436128 


87.25 


11200 


5255 


0.185 


0.219 


1, 2 


KO Aql 


AO + ... 


2.863965 


77.86 


9900 


4425 


0.144 


0.277 


3, 4 



1 Zola (1992), 2 Kreiner (2004), 3 Sahade (1945), 4 Mader & Angione (1996) 



sion features are presented in §5 and §6. Finally, the results 
are discussed in §7 and summarized in §8. 

2 OBSERVATIONS AND DATA REDUCTION 

The spectroscopic observations have been performed at the 
Catania Astrophysical Observatory with the echelle spec- 
trograph (FRESCO), which is fiber-linked to the 91-cm 
telescope. The data were obtained during two observing 
runs, with observing nights spread from 12 September to 
19 November 2003 and from 12 April to 29 July 2004. The 
spectra were obtained using the echelle configuration based 
on a 300 lines/mm echellette grating as cross-disperser ele- 
ment and an echelle grating with 79 lines/mm. A CCD cam- 
era with a thinned back-illuminated SITe chip of 1024 x 1024 
pixels (pixel size 24 fj,m) refrigerated to — 130° C by liquid 
nitrogen was used to acquire the data. The spectra have a 
resolution R ~ 21 000 and cover the wavelength range from 
4300 to 6650 A separated into 19 orders. We have obtained 
27 and 37 spectra for S Equ and KO Aql, respectively. The 
integration times for the variables were in the range 40-55 
minutes, and the S/N ratio reached was typically between 
50 and 80 at the continuum near the Ha, depending on the 
weather conditions. More information about the observa- 
tions can be found in Table 2. 

In addition to S Equ and KO Aql, we also observed some 
reference stars, whose spectral types are similar to those of 
the components of the variables. We have chosen the non- 
active and slowly-rotating standard stars 10 Tau (F9 IV-V) 
and a Boo (K1.5 III) to mimic the cooler components of 
S Equ and KO Aql, respectively. The bright, slowly-rotating 
stars, Vega (AO V, V r = -13.8 kms -1 ), /3 Vir (F9 V, V r = 
+4.6 kms -1 ), 5 Ser (F8 IH-IV, V r = +53.5 kms -1 ) and 
a Boo (K1.5 III, V r = —5.2 kms -1 ) have been chosen as 
templates for the radial velocity measurements of the hotter 



Table 2. Information about the spectroscopic observations 





Observing 


Number 


Exposure 


Number 


System 


Runs* 


of nights 


time (s) 


of spectra 


S Equ 


1 


8 


2400-3300 


11 




2 


15 


2400-3300 


16 


KO Aql 


1 


10 


2400-3300 


11 




2 


20 


2400-3300 


26 


Standard Stars 


Vega (A0V) 


1 


8 


15-30 


8 




2 


19 


15-30 


20 


10 Tau (F9IV-V) 


1 


2 


600 


2 


f3 Vir (F9V) 


2 


10 


300-600 


10 


5 Ser (F8III-IV) 


2 


2 


720 


2 


a Boo (K1.5III) 


2 


30 


20-30 


35 



*Run 1: 12 September - 19 November 2003 
Run 2: 12 April - 29 July 2004 



and cooler components and have been observed several times 
during each run. 

The data reduction was performed by using the 
ECHELLE task of IRAlf] package following the standard 
steps: background subtraction, division by a flat field spec- 
trum given by a halogen lamp, wavelength calibration using 
the emission lines of a Thorium- Argon lamp, and normal- 
ization through a polynomial fit to the continuum. 

The red part of the spectra is affected by many water 
vapor and O2 telluric lines. We have removed these lines 
at the Ha wavelengths using the spectra of Altair (A7 V, 
vsin i = 245 kms -1 ) and a Leo (B7 V, v sin i = 353 km s -1 ) 
acquired during the observing runs. These spectra have been 
normalized, also inside the very broad Ha profile, to provide 
valuable templates for the water vapor lines. Applying an 
interactive method that is described by Frasca et al. (2000), 

1 IRAF is distributed by the National Optical Astronomy, which is 
operated by the Association of University for Research in Astron- 
omy, inc. (AURA) under cooperative agreement with the National 
Science Foundation. 



the telluric lines were cleaned from the spectra of S Equ and 
KO Aql. 

3 RADIAL VELOCITIES AND ORBIT 
SOLUTIONS 

The cross-correlation technique (CCT) has been widely and 
successfully used, after the advent of electronic detectors, for 
the measurements of radial velocity (RV) reaching precisions 
as small as a few tens of ms -1 in particular cases, leading to 
the detection of giant planets around solar-type stars (e.g., 
Mayor & Queloz, 1995). 

By means of CCT, RVs of many binary systems were 
obtained with better precision and, in several single-lined 
binaries, the signal of the previously unseen companion has 
been detected. In the case of Algol-type binaries, the RV 
measurement of the cooler components is normally ham- 
pered by the strong luminosity contrast between the primary 
and secondary stars. The CCT has been proven to be very 
powerful for detecting the RVs of the secondary components 
of Algol systems (e.g., Holmgren 1988, Hill 1993, Khalesseh 
& Hill 1992, Maxted & Hilditch 1995). 

In the present work, RVs of both components of S Equ 
and KO Aql were derived by means of the CCT using the 
I RAF task FXCOR (e.g. Tonry & Davis 1979, Popper & 
Jeong 1994). The wavelength ranges were selected to exclude 
Balmer and Na I D2 lines, which strongly broaden the cross- 
correlation function (CCF) with their wings and can also 
be contaminated by the mass transfer and chromospheric 
activity. The spectral regions heavily affected by telluric ab- 
sorption lines (e.g. the O2 band A 6276-6315) were not used. 
The RVs of the components of S Equ and KO Aql, listed 
in Tables 3 and 4 together with their standard errors, are 
the weighted averages of the values obtained from the cross- 
correlation of each order of the target spectra with the cor- 
responding order of the standard star spectrum. The weight 
Wi = 1/crf has been given to each measurement. The stan- 
dard errors of the weighted means have been calculated on 
the basis of the errors (ai) in the RV values for each order 
according to the usual formula (e.g., Topping, 1972). The 
CTi values are computed by FXCOR according to the fitted 
peak height as described by Tonry & Davis (1979). 
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Spectra of the RV standard stars a Boo (K1.5 III) and 
(5 Vir (F9 V) have been used as templates for deriving the 
radial velocity of the secondaries of KO Aql and S Equ, 
respectively. An average high-S/N spectrum of Vega was in- 
stead used as template for the measurement of the RV of 
the primaries. As pointed out by Griffin (1999), it is not 
easy to achieve very accurate radial velocity measurements 
(accuracy better than 100 ms -1 ) for early type stars, due 
to the small number of suitable lines spread in a wide wave- 
length range and to other factors affecting the early-type 
line spectra. She found that real stellar spectra can be used 
as RV templates, because they are free from uncertainties in 
the tabulated wavelengths and have line shapes and inten- 
sities that are actually encountered. From cross-correlation 
of several high-resolution spectra of Vega with an average 
template built up with Sirius spectra, she found RV errors 
of about 250-350 ms -1 . These errors, possibly systematic, 
could be related to the known low-amplitude S Scuti-type 
variability of Vega (e.g., Samus et al. 2004) but are much 
smaller than the standard errors on the individual RVs of 
the primary components of S Equ and KO Aql (typically 
from 1 to 7 kms -1 ), and so they can be neglected. 

3.1 S Equ 

The RV curve of the primary component of S Equ has been 
obtained by Plavec (1966) based on plate spectra, which 
were taken by Petrie and Thackeray at Dominion Astro- 
physical Observatory and Radcliffe Observatory. The orbital 
solution obtained by them and, subsequently, by Lucy & 
Sweeney (1971) on the same very scattered data is charac- 
terized by a relatively large eccentricity (e « 0.15). This 
rather high value is not supported by the light curves with 
equally spaced minima of the same duration observed by 
Catalano & Rodono (1968) and Zola (1992). Our RV data 
definitely exclude an eccentricity larger than about 0.02, and 
so a circular orbit was assumed in the following analysis. Up 
to now, there has been no published RV curve of the sec- 
ondary component, to our knowledge. 

In Figure 1, the measured RVs and the associated error 
bars of the components of S Equ are plotted as a function 
of the orbital phase. Filled and open circles represent the 
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Table 3. Radial velocity measurements of S Equ 



HJD 


Orbital 


Vh 






Vc 




2 450 000+ 


Phase 


(km s~ 




(km s - 


) 


2895 4598 





1907 


— 81. 7± 


2 


6 


122 4+1 5 


5 







4781 


— 55.8+ 


3 









291 4581 





5556 


— 45. 5± 


6 


9 






9Q99 Ifilf) 





0202 


— 79.8+ 


5 


6 


—32.3+ 3 


X 


2922 4364 





0416 


— 66.1+ 


6 


5 


— 12.6+ 4 


2 


9Q93 40^8 





3237 


— 72.8+ 


9 


1 


121.0+ 7 




2924 41 77 





6182 


— 37. 3± 


A 


9 


— 198.9+ 4 


7 







2207 


— 79.3+ 


7 


1 


131.6+ 4 


7 


9Q^4 2877 





4906 


— 50. 2± 


5 


1 


—39.1+15 


o 


9Q34 334^ 





5042 


— 49.9+ 


3 


8 


—44.6+ 10 









5122 


— 49. 3± 


1 


6 


— 102 8+ 8 


x 


31 ^ 





8779 


— 38. 2± 


1 


3 


— 184.2+ 2 


9 


31 57 5783 





4738 


— 56. 1± 


2 














7591 


— 28. 5± 


1 


9 


—247.2+ 3 


4 


3160.5341 





3340 


-74.5± 


2 


6 


116.7+ 3 


5 


3166.5869 





0955 


-68. 5± 


1 


9 


55.7+ 3 


1 


3167.5542 





3770 


-70.9± 


1 


4 


72.0+ 2 


6 


3168.5489 





6665 


-32.2± 


2 


3 


-227.4+ 3 


7 


3170.5519 





2494 


-79.3+ 


1 


7 


135.7+ 3 


7 


3171.5905 





5517 


-49. 5± 


1 


9 


-123.7+ 8 





3185.5247 





6069 


-37.2± 


2 


5 


-185.1+ 4 


5 


3186.4374 





8725 


-36. 7± 


1 


7 


-184.5+ 3 


8 


3202.5972 





5754 


-43. 2+ 


1 


5 


-152.0+ 5 


2 


3203.4302 





8178 


-30.9± 


1 


6 


-226.0+ 4 





3213.4549 





7353 


-29.6± 


1 


9 


-243.5+ 5 





3214.4695 





0290 


-65.0± 


1 


fi 


-21.6+ 2 


9 


3214.5830 





0636 


-63. 7+ 


1 


7 


22.4+ 2 


2 



primary and secondary components, respectively. The or- 
bital solution (dashed line for the cooler and solid line for 
the hotter component) was determined, assuming circular 
orbits, by a least-squares fit which reproduces very well all 
the observed RVs. We report in Table 5 the orbital param- 
eters of the system derived from the solution of our radial 
velocity curves and compare them with with the values de- 
rived by Plavec. 

In the orbital solution of Plavec (1966) given in Table 
5, the systemic velocity, 7, and the semi-amplitude of the 
RV curve, Ki, were found to be —47.95 ± 0.6 kms -1 and 
23.4 + 0.8 kms -1 , respectively. 



3.2 KO Aql 

Sahade (1945) obtained the first RV measurements and or- 
bital solution for the primary component of KO Aql based 
on low-resolution spectra. He was not able to detect the 
secondary component due to its very low contribution to 
the spectrum of the system observed at blue wavelengths. 
Then, Lucy & Sweeney (1971) performed a new orbital so- 



S Equ 

, 1 , , , , 1 , , , , r 




, , , I , , , , I , , , , I , , , , I , , , , I , , , 1 

. 0.25 0.50 0.75 1.00 

Phase 

Figure 1. Radial velocity curves of a double-lined orbit for S Equ 
plotted as a function of orbital phase. Filled and open circles 
represent the RVs of primary and secondary components, respec- 
tively. Solid line represents the orbital solution for the more mas- 
sive component, while the dashed line for the less massive com- 
ponent. 

lution for the primary component with results similar to 
those of Sahade (1945). We have obtained high-resolution 
spectra that allowed us to determine the orbital parameters 
of the cooler component for the first time and those of the 
hotter one with a better precision than in previous stud- 
ies. The RVs were measured from the observed spectra us- 
ing the cross-correlation technique as described above. Two 
template spectra, Vega for primary component and a Boo 
for secondary component, have been used for the CCT. The 
final RVs for KO Aql are listed in Table 4. The larger er- 
rors in the velocities of the secondary (up to 3-4 times the 
primary's ones) are due to the very low contribution of the 
secondary star to the observed spectrum and to its high ro- 
tational velocity that broaden and weaken its spectral lines. 
The fit of an eccentric solution to the observed RVs of both 
components provided very low, not-significant values of ec- 
centricity (e w 0.01). Thus, a circular orbit was assumed for 
the fit of the observed RV curves, obtaining the elements 
reported in Table 5, where the ephemeris was adopted from 
Kreiner (2004). 

A plot of the observed RVs together with the fitted RV 
curves is shown in Fig. 2. The orbital solutions determined 
by Sahade (1945) for only the primary component and that 



Table 4. Radial velocity measurements of KO Aql 
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KO Aql 



HJD 


Orbital 








V c 




2 450 000+ 


Phase 


(km s 




) 


(km s 




2895 4293 





4696 


— 7.1+ 


4 


2 






98Qfl 111D 
jl. O iJ u . >_> >_> o u 





7851 


42.7+ 


3 


4 






291 3800 





6897 


39. 2± 


2 


x 


— X64.4+ 


7.2 


9Q99 971 4 





8416 


33. 5± 




3 


— 142.3+ 


5.8 


2923 2988 





2003 


— 32.3+ 


x 


2 


X65.3+ 


3.6 


9Q94 1989 





5598 


12.8+ 


X 


x 


—62.3+ 


7.8 


2925 3135 





9038 


22.4+ 


x 


5 


—85.1+ 


4.7 


9Q11 1917 





7006 


39.5+ 


x 


5 


— 165.4+ 


6.1 


2934 2506 





0242 


— 15.2+ 


x 


5 


21.0+ 


3.8 


9Q^Q 9777 





7625 


41.9+ 


3 


g 


— 173.8+ 


9.1 


3125 5480 





8164 


37.6+ 


x 


5 


— 155.4+ 


3.6 


11 9fi dQR^ 





1483 


—29.0+ 


x 


3 


146.6+ 


2.7 


3127 5205 





5051 


0.9+ 


x 


x 






11 17 ^019 





9906 


8.0+ 


2 


2 


—7.5+ 


2.1 


31 37 5437 





0047 


1.2+ 


2 


x 


—0.9+ 


2.8 


iJIOO.'JOUU 





3493 


—26.5+ 


3 


4 


145.3+ 


1.7 


•! 1 30 ^fi25 





7096 


36.5+ 


x 


4 


— 170.9+ 


4.3 


11 40 489^ 





0308 


— 12.4+ 


x 


4 


31.9+ 


3.9 


3140 5fifi9 





0603 


— 11.0+ 


x 


4 


66.8+ 


5.7 


11 41 4Q7Q 





0837 


— 17.0+ 


x 


5 


86.8+ 


4.0 


3143 c;831 





1134 


—23.7+ 


x 


5 


109.7+ 


2.2 


11 ^9 4QD c i 





2235 


—36.5+ 


x 


x 


173.4+ 


5.4 


31 cjo 5880 





2575 


—34.6+ 


x 


x 


181.0+ 


3.1 


3155.5687 





2892 


-34.7+ 





9 


166.5+ 


4.9 


3156.5837 





6526 


32.9+ 


1 





-144.5+ 


3.0 


3158.5267 





3310 


-31.7+ 





9 


154.4+ 


5.4 


3161.5530 





3877 


-20.6+ 


1 


5 


123.8+ 


3.3 


3166.5406 





1291 


-25.0+ 


1 


2 


124.2+ 


6.0 


3167.4317 





4403 


-9.7+ 


1 


9 


78.0+ 


6.2 


3167.5478 





4808 


-0.1+ 





7 


26.9+ 


2.7 


3168.4281 





7882 


39.9+ 


1 


3 


-165.1+ 


3.3 


3168.5833 





8423 


31.7+ 





9 


-141.0+ 


3.1 


3170.5860 





5416 


10.2+ 





7 


-51.2+ 


3.7 


3171.5641 





8831 


27.8+ 





9 
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2.2 


3216.5501 





5901 


22.4+ 


1 





-98.4+ 


2.7 



obtained in this work (see Table 5 for comparison) agree 
very well with each other. 



4 ROTATIONAL VELOCITIES 

The synchronization time is shorter than the duration of 
the semi-detached phase for Algol-type binaries by two or 
three orders of magnitude. In this case, one should expect 
that the number of Algol systems with a non-synchronous 
rotation must be much smaller than that of synchronous 
systems (Glazunova 1999). However, the measured rotation 
rates of the primary components of Algol systems indicate 
that most of them are rotating faster than in the case of 
synchronization (Van Hamme & Wilson 1990). The mass 
transfer process from the secondary components to the pri- 




Figure 2. Radial velocity curves of the components of KO Aql 
plotted as a function of orbital phase. The symbols and lines have 
the same meaning as in Fig. 1. 



mary ones can cause a strong increases in the rotational 
velocity of the mass-gaining components because the mat- 
ter flowing through the Li Lagrangian point also transfers 
angular momentum (see, e.g., Giuricin et al. 1984, De Greve 
at al. 2006, Van Resbergen et al. 2006). 

Van Resbergen et al. (2006) have stated that a "conser- 
vative" evolution (without loss of mass and angular momen- 
tum) fails to explain the observed mass-ratios (q = 0.4-1) 
while models with some mass and angular momentum loss 
from these systems can better explain the observed mass- 
ratios. They have proposed a scenario in which the gainer is 
spun up by the incoming mass from the donor. The gainer 
spins up - sometimes up to the critical velocity - and spins 
also down due to tidal interaction. Even when the gainer 
has the critical velocity this kinetic energy is not sufficient 
to drive mass out of the system, as a consequence of the 
virial theorem. The extra-energy needed is given by the 
accretion luminosity which is created in a hot spot of the 
gainer's surface. As a consequence, they found that, in mas- 
sive Algol systems (with a donor starting with an initial 
mass ^ 5Mq), the combined action of rapid rotation and 
accretion luminosity can overrun the binding energy of ma- 
terial on the gainer's surface. The less-massive Algol systems 
(like, e.g., RZ Cas) are instead likely obtained through an 
evolution without loss of mass but with a considerable loss 
of angular momentum (e.g. due to magnetic braking). 
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Table 5. Spectroscopic orbital elements of S Equ and KO Aql. 







S Equ 


KO Aql 


Parameter 


Plavec (1966) 


This study 


Sahade (1945) 


This study 


To (HJD) 


2437568.345 


2452503. 086 1 


2426585.442 


2452501. 707 1 


P ori) (day) 


3.4360726 


3.436128 1 


2.863844 


2.864068 1 


e 


0.14 


0.01 ±0.01 2 


0.02 


O.OliO.01 2 


7 (km s _1 ) 


-47.95 


-54.5± 0.5 


-2.7 


1.9 ± 0.1 


Ki (km s" 1 ) 


23.4 


25.1 ± 0.6 


37.8 


37.8 ± 0.2 


K 2 (km s" 1 ) 




192.1 ± 1.2 




175.1 ± 0.9 


ai sin i (10 6 km) 


1.06 


1.19 ± 0.03 


1.5 


1.49 ± 0.01 


a2 sin i (10 6 km) 




9.08 ± 0.06 




6.90 ± 0.03 


Mi sin 3 i (Mq) 




3.23 ± 0.06 




2.36 ± 0.03 


M2 sin 3 i (Mq) 




0.42 ± 0.01 




0.51 ± 0.01 


q (=M 2 /Mi) 




0.131 ± 0.003 




0.216 ± 0.002 



The ephemeris was adopted from Kreiner (2004). 
2 e = was assumed. 



Therefore, the v sin i of the accreting components repre- 
sents observational evidence for mass transfer and accretion 
phenomena. Moreover, the values of v sin i are needed to 
produce correct synthetic spectra. 

We measured the projected rotational velocity (v sin i) 
of the "mass gainers" (i.e. hotter components) of S Equ and 
KO Aql by using the CCF analysis. For this purpose, we 
obtained a calibration relation between the full width at 
half maximum (FWHM) of the CCF peak and the v sin i 
with artificially broadened spectra of slowly-rotating B9-A0 
stars (Vega, 21 Peg, and HR 6096) that were acquired in 
the same observing runs and with the same instrumental 
setup of S Equ and KO Aql observations. Then, the rota- 
tion velocities of the hotter components of the systems were 
determined by converting the FWHM of its CCF peaks into 
■usini through the aforementioned calibration. We found, as 
weighted averages of all measurements, fsini = 52.4 ±4.4 
kms^ 1 and 41.0 ±2.0 kms" 1 for S Equ and KO Aql, re- 
spectively. The rotational velocity for S Equ is in agree- 
ment, within the errors, with the value of wsini = 47.4 ± 1.4 
kms -1 determined by Mukherjee et al. (1996). 

We have also calculated the rotational velocities of both 
components of the systems in the hypothesis of synchronous 
rotation, finding '■-'sync Sill % = 40km s _1 and 30 kms -1 for 
the hotter components of S Equ and KO Aql, respectively 
(see, Table 9). The measured usini values for the gainers 



indicate that they are spinning faster than the synchronous 
rate. 

5 SYNTHETIC SPECTRA AND Ha 
DIFFERENCE PROFILES 

The emission and absorption features produced by accretion 
around the hotter component in short-period Algol systems 
are generally weak relative to the continuum flux of the spec- 
tra. In order to emphasize the effect of mass transfer (and/or 
other physical processes such as magnetic activity in cool 
components) in Ha, it is necessary to remove the strong 
photospheric contribution of both components from the ob- 
served spectra. Hence, we have produced synthetic spec- 
tra for the hotter components built up with Kurucz (1979) 
model atmospheres. Standard stars' spectra have been used 
for the cooler components. For this purpose, the low-activity 
stars 10 Tau (F9 IV- V) and a Boo (K1.5 III) for S Equ and 
KO Aql, respectively, were observed in the same observing 
runs and with the same instrumental setup. 

Before forming the difference profiles (observed - syn- 
thetic) for the Ha line, we have checked the physical pa- 
rameters (T e ft and log g) of the hotter components reported 
in the literature. This task was accomplished by fitting Ku- 
rucz (1979) model atmosphere spectra to the observed ones 
in spectral regions where hot components overwhelm the 
secondary ones and mass transfer effects are not visible. We 



S Equ 




0.6 



4520 4530 4540 4550 4560 4570 4580 4590 
Wavelength (A) 
KO Aql 




0.6 g = 

4520 4530 4540 4550 4560 4570 4580 4590 
Wavelength (A) 

Figure 3. Comparison between the observed (thin line) and the 
synthetic spectra (thick line), normalized to a continuum of unity 
in the blue spectral region that gives the better match for S Equ 
(a) and KO Aql (b). 

have chosen the blue region around 4500 A and H7. We 
have also taken into account the small contribution (from 
3 to 7 %) of the cooler components at these wavelengths 
by adding observed spectra of cool standard stars properly 
weighted, broadened and Doppler-shifted. The temperature 
of the primary component of S Equ found using this method, 
Ti = 11, 500 K, is only slightly higher (300 K) than that 
given by Zola (1992). For KO Aql, we find the same tem- 
perature (Ti = 9900 K) found by Mader & Angione (1996). 

The stellar parameters resulting from these fits are 
listed in Table 6, which contains also other parameters used 
to create the synthetic spectra for S Equ and KO Aql. Some 
examples of observed and synthetic spectra are displayed 
in Fig. 3. Then, we have formed the composite synthetic 
spectra in the Ha region to be compared with the observed 
ones and to produce difference profiles. Effective temper- 
atures and absolute radii (in Table 9) were used to calcu- 
late the contribution of primary and secondary components, 
whose spectra were also rotationally broadened and wave- 
length Doppler-shifted, before being summed. The synthetic 
profiles were then subtracted from the observed ones, and 
the resulting difference profiles were examined. 
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6 Ha EXTRA-EMISSION AND ABSORPTION 

The Ha residual profiles of Algol systems generally display 
extra absorption and/or emission features commonly inter- 
preted as the result of Roche-lobe overflow, accretion struc- 
tures around the hotter components, and/or due to the chro- 
mospheric activity of the cooler components. In short-period 
systems the contribution of the active chromosphere of the 
rapidly-rotating cooler components can produce significant 
effects on the Ha line profile (Richards & Albright 1999). 
The two Algol systems with relatively short periods and 
with late-type secondary components that we are study- 
ing present both effects, namely mass transfer and chromo- 
spheric activity. 

6.1 S Equ 

The Ha difference profiles of S Equ display extra absorption 
and emission that are variable with orbital phase. These 
spectral features originate from accretion structure around 
the primary, flowing gas in the region between the two stars 
and by the impact of the gas stream on the primary compo- 
nent (Richards & Albright 1999, Vesper et al. 2001, Richards 
2004). Ha emission closely associated with the late- type sec- 
ondary component and attributed to chromospheric activity 
has been also detected by Richards (2001) in the Ha tomo- 
gram of S Equ. 

In this work, we have analyzed Ha difference profiles 
by selecting the spectra with higher S/N ratio. Examples of 
the observed, synthetic and difference Ha profiles of S Equ 
are shown in Fig. 4. Several features could be readily noted. 
Some extra emission is present at nearly all phases except 
near primary eclipse. These emission components with single 
or double peaks are more clearly visible around the quadra- 
tures and the stronger peaks, which are red shifted, can be 
seen near phase 0.25. Although the shape of the features at 
phase ~ 0.75 is similar, the strength is slightly weaker. The 
extra-absorption is more prominent just before and after pri- 
mary eclipse and it is still present at other phases though 
with a decreasing strength. 

Considering all the residual profiles obtained at differ- 
ent orbital phases, a general structure, consisting of a weak 
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Table 6. The model parameters and standard stars for S Equ and KO Aql to create synthetic spectra. 





Hot Component 


Cool Component 


System 


T log g vsini 


Standard 


v sin i 




(K) (kms" 1 ) 




(kms" 1 ) 


S Equ 


11500 4.1 52 


10 Tau (F9 IV) 


■18 


KO Aql 


9900 4.3 41 


a Boo (K1.5 III) 


58 



27/07/2004 



09/06/2004 



13/06/2004 



03/06/2004 



15/07/2004 



01/06/2004 



29/05/2004 



Equ - Ha 

0.2 

o. 
o. 

-o. 

0. 

o. 

0. 
-0. 

0. 

0. 

0. 
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0. 
0. 

0. 
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0. 

0. 

0. 
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0. 

0. 

0. 
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0. 
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0. 

0. 

0. 
-0. 
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Figure 4. Left. Sample of observed (thin line) and the synthetic 
spectra (thick line) for S Equ at different orbital phases. The 
date of observation is marked in each panel. Right. The difference 
spectra showing H« extra-absorption and emission features. The 
orbital phase is reported in the upper right corner of each panel. 
The vertical dashed lines in the lefts and right boxes represent 
the laboratory wavelength (Aq = 6562.82) of Ha line. 



absorption (A) and broad, single- or double-peak emission 
(B and C), is met (see Fig. 5). 

With the aim of determining the sources of such fea- 
tures, we measured the equivalent widths (EWs) and radial 



velocities (RVs) for the absorption and emission components 
in the difference profiles by means of multiple Gaussian fits. 
These values are reported in Table 7. Moreover, in some 
spectra we have detected other features which seem to be 
unrelated to the ones already identified as A, B, and C, on 
the basis of their RV. These unidentified features are referred 
to simply as "absorption" or "emission" in Table 7. 

The EW and RV values for the A, B, and C features 
are plotted versus orbital phase in Fig. 6 with asterisks, plus 
signs, and triangles, respectively. The RV curves of the pri- 
mary and secondary stars are also over-plotted in the same 
figure with full and dotted lines, respectively. The sinusoidal 
fit of the RV variation of the B emission feature is displayed 
by a dashed line. 

Extra Absorption Feature (A). The visibility and ve- 
locity information, which can be deduced from EW and RV 
variations, provide constraints about the location of sources. 
Ha extra absorption is clearly visible just before and after 
primary eclipse (Fig. 4). It can be also easily followed in 
Fig. 5, which shows the EWs and RVs of the A-absorption 
feature. The EWs indicate that the absorption is absent at 
phases far from the primary eclipse. The RVs of this extra- 
absorption feature closely follow the RV curve of the hotter 
component suggesting that it may be connected with the 
mass-gainer component (Fig. 6a). This feature could be the 
result of circumstellar matter arising from the interaction 
between the gas stream and the environment around the 
hotter star. The EW variations also indicate that the mat- 
ter should be denser around the impact point, where the gas 
stream interacts with the primary star, which is more visi- 
ble around primary eclipse. The strong extra-absorption ob- 
served just before the primary minimum could also be the ef- 
fect of the gas stream itself that is projected against the disk 
of the hotter component at these phases. However, the extra- 
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Table 7. The EWs and RVs of absorption and emission features 
in Ha difference profiles for S Equ. 



HJD Orbital Type 1 BW V r 

24 50000+ Phase (A) (km s -1 ) 



-61 
408 
282 
163 

20 
125 
-43 
180 
-30 
121 
244 
-98 

16 
121 
-73 
1 

74 
-23 
1 10 
225 
-87 

36 
137 
123 
165 
-42 
215 
-35 
158 

93 
-89 
199 
-72 



1 The features labeled as "absorption" and "emission" have different RV 
variations from those listed as A, B and C. 
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Figure 5. A difference Ho profile for S Equ which displays all 
the prominent features; the extra-absorption (A) the broad emis- 
sion component (B) and the narrower emission connected to the 
secondary star (C). 



absorption observed just after primary minimum and its RV 
variation make the impact region the more likely source of 
this feature. This would also explain the light deficit ob- 
served by Catalano & Rodond (1968) near 0.9 phase. 

Extra Emission Feature (B). The EWs show that the 
emission profile (B) is not seen during the eclipses, while it 



is more prominent around the quadratures (see Figs. 4 and 
6b). Therefore, it may be related to the circumstellar mat- 
ter located in the space between the primary and secondary 
star, that is occulted at both eclipses (see Fig. 7). As shown 
in Fig. 6a, the RVs of the B emission feature follow the or- 
bital motion of the cooler component without any significant 
phase shift but with a smaller amplitude. In the hypothesis 
that the emitting region is at rest in the rotating reference 
frame of the system, we can deduce that it lies between 
the center-of-mass of the system (CM) and the secondary 
star and we can also evaluate its location, analogously to 
what done for HR 7428 by Marino et al. (2001). The pro- 
jected distance of the region from the CM can be found as 
o.b = cl c {Kb / K c ), where a c = a/(l + q) is the semi-major 
axis of the cool star orbit and q and a are the mass ratio 
and system separation, respectively. Here, K c and K B are 
the semi-amplitudes of the RV variations for the cool star 
and for the emitting region, respectively. Using the orbital 
solution of S Equ given in Table 5 and the semi-amplitude 
of the RV variation of the B emission feature estimated as 
Kb — 71 kms -1 , we find that the emitting region should be 
located at about 3 x 10 6 km from the barycenter toward the 
cool component. Thus, feature B could be part of the gas 
stream (see Fig. 7). 

Extra Emission Feature (C). Although the emission pro- 
file (C) is visible at nearly the same orbital phases where the 
B-type emission is observed, its narrower width (Fig. 5) and 
its RV variation (Fig. 6a) clearly indicate that its origin is 
different. The velocity of the C feature follows the RV curve 
of the secondary companion. Therefore, the emission seems 
to be closely related to the cool component and it could be 
produced by its chromosphere. 

A schematic representation of the S Equ system with 
the Roche-lobe geometry, the gas stream, and the extra- 
emission/absorption regions is shown in Fig. 7. 

6.2 KO Aql 

We have subtracted the synthetic spectrum from each 
KO Aql spectrum using the method described in the pre- 
vious section in order to produce residual Ha profiles to be 
analyzed. A sample of Ha profiles is shown in Fig. 8. In the 
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Figure 6. The variations of RVs (upper panel) and EWs (lower 
panel) of extra-absorption and extra-emission features observed 
in H« for S Equ. The continuous and dotted lines are the RV 
solutions for the primary and secondary component, respectively. 
The dashed sinusoid represents the RV variation of the broad 
emission region (B) with a semi-amplitude Kb — 71kms _1 . 



SEqu 




Figure 7. Schematic representation of the S Equ system 
with the Roche-lobe geometry, the gas stream, and the extra- 
absorption/emission regions (A, B, and C). The barycenter, which 
is inside the hotter component, is marked by a cross. 



residual profiles is apparent an absorption feature (A), which 
is stronger around the primary eclipse, and an emission pro- 
file (B) clearly seen around the two quadratures. We have 
measured EWs and RVs of these features to get information 
about their origin, as we did for S Equ. These values are 
listed in Table 8 and are plotted against the orbital phase 
in Fig. 9. However, features with RV not consistent with A 
and B are visible in a few spectra. The EWs and RVs for 
these features are not plotted in Fig. 9, but their values are 
reported in Table 8 and simply labeled as "absorption" or 
"emission" features. 

Although the RV of the extra-absorption (A) seems to 
follow the RV curve of the mass-gainer (hotter) component 
at some phases (such as between 0.25 and 0.35), it dis- 
plays a much larger amplitude of RV variation compared 
to the primary star, especially just before and after pri- 
mary eclipse (see Fig. 9). The EWs indicate that the ab- 
sorption feature becomes stronger from phase 0.7 to 0.9 and 
is visible until phase 0.2. The maximum value of equivalent 
width, EW»;0.6A, is observed just before primary eclipse, 
at about 0.9 phase. This suggests that the absorption is 
likely produced by circumstellar matter around the primary 



component or to the gas stream itself that is projected onto 
the hotter star and would give rise to a red-shifted absorp- 
tion at these phases, as it is observed. However, this can- 
not explain the blue-shifted absorption observed at about 
0.1 phase. Anyway, the variability and the visibility of the 
absorption feature (A) indicate that the source of this extra- 
absorption may be connected with the stream-star impact 
region. 

To estimate the location and the origin of the extra- 
emission feature (B), we examine the variation of its RVs and 
EWs in Fig. 9. The EWs have a value of about 0.2 A and no 
relevant change with orbital phase. The RV of this emission 
feature closely follows that of the secondary star, indicating 
a possible chromospheric origin. Unlike S Equ, for KO Aql 
there is no clear evidence of a broad extra-emission feature 
attributable to the mass transfer phenomenon. A schematic 
representation of the KO Aql system with the Roche-lobe ge- 
ometry, the gas stream, and the extra-emission/absorption 
regions is shown in Fig. 10. 
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Table 8. The EWs and RVs of absorption and emission features 
in Ha difference profiles for KO Aql. 
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Figure 8. Left. Sample of observed (thin line) and the synthetic 
spectra (thick line) for KO Aql at different orbital phases. The 
date of observation is marked in each panel. Right. The difference 
spectra showing Ha extra-absorption and emission features. The 
orbital phase is reported in the upper right corner of each panel. 
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29 


3166.5406 


0.1291 


A 


0.289 


-109 


3166.5406 


0.1291 


B 


-0.256 


145 


3166.5406 


0.1291 


Absorption 


0.096 


-88 


3168.4281 


0.7882 


A 


0.223 


125 


3168.4281 


0.7882 


B 


-0.304 


-152 


3168.5833 


0.8423 


A 


0.433 


144 


3168.5833 


0.8423 


B 


-0.137 


-142 


3171.5641 


0.8831 


A 


0.323 


151 


3171.5641 


0.8831 


B 


-0.120 


-111 


3216.5501 


0.5901 


A 


0.144 


170 


3216.5501 


0.5901 


Absorption 


0.063 


-189 



1 The features labeled as "absorption' 1 '' and "emission" have different RV 
variations from those listed as A, B and C. 



7 DISCUSSION 

According to the catalogue compiled by Budding et al. 
(2004), the number of Algol-type binaries with RV curves for 
both components and solved orbits is very limited. Accurate 
values of stellar and orbital parameters for as many bina- 
ries as possible are strongly needed for the study of several 
phenomena such as mass transfer and accretion, magnetic 
activity on late-type components, and angular momentum 
evolution. 

In this study, we have presented new radial velocity 
curves for the hotter primary components of S Equ and 
KO Aql with much better precision and phase coverage than 



previous works. Moreover, we have obtained the first RV 
curves for the late-type secondary components of both sys- 
tems. 

Thus, the first complete orbital solutions and mass de- 
termination have been made for both systems. Prom the 
photometric and spectroscopic elements given in Tables 1 
and 5, respectively, we derived the absolute parameters of 
S Equ and KO Aql, which are listed in Table 9. The masses 
are determined to better than 2%, and the radii with a 4% 
precision. 

From the effective temperatures and radii we have cal- 
culated the luminosities of the components of both systems 
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Figure 9. The variations of RVs and EWs of Ha difference profiles 
for KO Aql. As before, the continuous and dotted lines represent 
the RV solutions for the primary and secondary, respectively. 



KO Aql 




Figure 10. Schematic representation of the KO Aql system 
with the Roche-lobe geometry, the gas stream, and the extra- 
absorption/emission regions (A and B). The barycenter, close to 
the hotter component, is marked by a cross. 



and their absolu te bolome tric magnitudes, adopting for the 



Sun M®j =4.74 l|Coxll2000l ). The absolute magnitudes in the 
V band, Mv, were calculated from the bolometric magni- 
tudes by applying the bolometric corrections appropriate 
for the temperatures of our stars (Flower 1996). The com- 
bined absolute magnitude (including both components) is 



have also calculated the combined B — V color expected from 
the sources by means of the calibration T e g versus B — V 
provided by Flower (1996), finding (B - Vf ot = -(T032 for 
S Equ and (B - V)l ot = +Cr046 for KO Aql. 

From the absolute and apparent V magnitudes we can 
easily derive the distance modulus, if we know the reddening 
to the sources. 

For S Equ, we have used the outside-of-eclipse values 
V = 8 m 37, B — V = (F072 (ESA 1997), that lead to a color 
excess E(B - V) ~ m 104 and an extinction Ay ~ (T324, 



with R = A V /E(B - V) = 3.1 l|Savage fc Mathis! Il979l ) 



A distance of 390 ± 40 pc is found, in excellent agreement 
with the astrometric determination of 400 ± 210 pc from the 
Hipparcos mission (ESA, 1997). 

For KO Aql, with V = 8 m 42 and B — V = (T114 (ESA 
1997), we find E{B - V) ^ m 068 and Av ~ m 212. The 
distance derived in the same way is 240 ± 20 pc, in good 
agreement with the astrometric distance of 265 ± 75 pc from 
the Hipparcos mission (ESA, 1997). 

The positions of the primary and secondary components 
of S Equ and KO Aql on the Hertzsprung-Russel (HR) di- 
agram are shown in Fig. 11, together with the theoretical 
evolutionary tracks for the more massive components, as 
computed by Girardi et al. (2000), and with the conserva- 
tive evolutionary track for the mass-gainer star belonging to 
a system with a total mass of Mi + Mi = 4.2 Mq and an 
initial orbital period of 1.2 days calculated by De Loore & 
Van Rensbergen (2005^. 

Ibanoglu et al. (2006) have pointed out that semi- 
detached systems (SDBs) display relevant differences com- 
pared to the detached binaries (DBs). The primary com- 
ponents, mass-gainers of SDBs, seem to be normal main- 
sequence stars. The cooler ones, the losers, are instead stars 
evolved off the main sequence which are all oversized and 
over-luminous for their masses. A relevant difference be- 
tween the primary components of SDBs and DBs has been 
found by Ibanoglu et al. (2006) who determined the mass- 
luminosity relation (MLR). They have taken into account 
74 DBs and 61 SDBs with well-determined absolute param- 



M 



0™095 for S Equ and M^ x = l m 326 for KO Aql. We 2 |http://we.vub.ac.be/astrofys/binarytracks/index.htrnT1 
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Table 9. Absolute parameters and related quantities for S Equ and KO Aql. vs'mi is the observed projected rotational velocity for both 
components. v syn sini is instead the projected rotational velocity of the primary components calculated assuming synchronous rotation. 





S Equ 


KO Aql 


Parameters 


Primary 


Secondary 


Primary 


Secondary 


Mass (Mq) 


3.24 ± 0.03 


0.42 ± 0.01 


2.53 ± 0.05 


0.55 ± 0.01 


Radius (Rq) 


2.74 ± 0.09 


3.24 ± 0.10 


1.74 ± 0.07 


3.34 ± 0.07 


log T (K) 


4.049 a 


3.72 a 


3.996 b 


3.646 6 


log L (Lq) 


2.02 ± 0.10 


0.86 ± 0.10 


1.41 ± 0.06 


0.56 ± 0.06 


log g (cgs) 


4.07 ± 0.01 


3.04 ± 0.01 


4.36 ± 0.03 


3.13 ± 0.01 


a(R ) 


14.8 


± 0.1 


12.05 


± 0.05 


v ayn sini (km/s) 


40.3 


47.7 


30.1 


57.7 


v sini (km/s) 


52.0 ± 4.0 




41.0 ± 2.0 




M boi (mag) 


-0.32± 0.11 


2.60 ± 0.15 


1.20 ± 0.13 


3.28 ± 0.15 


M, (mag) 


0.19 ± 0.11 


2.81 ± 0.15 


1.43 ± 0.13 


3.94 ± 0.15 


Distance (pc) 


390 


± 40 


240 


± 20 




a Zola (1992), 


^Mader & Angione 


(1996). 





eters and derived for the primaries of SDBs a MLR rela- 
tion Li oc Mf' 20 , while the luminosities of the primaries of 
DBs are correlated with their masses as M?- 92 . These rela- 
tions show that the hotter components of SDBs are under- 
luminous relative to the main-sequence stars of the same- 
mass. 

This trend is also followed by the hotter components 
of S Equ and KO Aql which are clearly underluminous 
compared to normal main-sequence stars of the same mass 
(Fig. 11). This behaviour is not dependent on the theoretical 
tracks, because we find exactly the same results by adopt- 
ing the evolutionary tracks of Pols et al. (1998), and so the 
trend could be due to the mass that is being acquired by 
these stars by Roche-lobe overflow from the donor. 

Indeed, the evolutionary track of the mass-gaining com- 
ponent for an Algol system with Mi + M 2 = 4.2 Mq (dot- 
ted line in Fig. 11) shows that, if the mass-transfer is still 
ongoing, the position of the hotter component can be signif- 
icantly different from that of a normal main-sequence star 
with the same mass. Unfortunately, De Loore & Van Rens- 
bergen (2005) did not calculate evolutionary tracks for sys- 
tems with the same total mass as S Equ (3.66 Mq) and 
KO Aql (3.03 Mq), but this result indicates that the lumi- 
nosity discrepancy could be explained by an evolutionary 
effect. 

S Equ and KO Aql are members of the group of short- 
period Algols (P or 6 < 5-6 days), in which the mass trans- 




KO Aql B 



/AM.": 



0.0 L 



4.1 4.0 3.9 3.8 3.7 3.6 

Figure 11. HR diagram of the primary (A) and secondary (B) 
components of S Equ (dots) and KO Aql (diamonds). The ZAMS 
and the post-main-sequence evolutionary tracks (full lines) for the 
masses of the primary components of S Equ and KO Aql from 
Girardi et al. (2000) are also shown for comparison. The evolu- 
tionary track of the mass-gainer component of a semi-detached 
system with a total mass of 4.2 Mq according to De Loore & Van 
Rensbergen (2005) is displayed with a dotted line. The thick grey 
line marks the portion of this track for which the mass of the 
gainer is comprised between 3.1 and 3.4 Mq. 



fer and accretion process are more complicated than in the 
long-period Algols (P or b > 5-6 days) because the gas stream 
from the donor to the gainer directly strikes the star. This 
stream-star interaction results in a variety of accretion fea- 
tures. Generally, an accretion structure in these systems, 
which is called an accretion annulus (Richards et al. 1996), 
can be produced. This structure is an asymmetric, some- 
times clumpy distribution of gas around the primary compo- 
nent (Richards 2001). The observed properties of the mass 
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Figure 12. The r-q diagram for classical Algols, where r is the 
fractional radius of the primary star and q is the mass ratio of 
the system. The curve of u> m i n represents the minimum distance 
between the gas stream and the center of the primary star, while 
the curve of u)d corresponds to the distance of closest approach of 
the stream. Notice that S Equ and KO Aql lie between the two 
curves, indicating unstable and transient accretion structures. 

transfer and circumstellar material in short-period Algols 
have been mostly found from spectroscopic studies, espe- 
cially from variations of Balmer-line profiles. The identified 
regions include the gas stream (Richards et al. 1995), a tran- 
sient accretion disk (Kaitchuck & Honeycutt 1982), and a 
denser region near the impact site called a localized region 
(Richards 1992, Richards et al. 2000). 

The type of accretion structures can be predicted by 
using the ballistic calculations of Lubow & Shu (1975) and 
the ri — q diagram, a plot of radius of the mass gainer versus 
the mass ratio (q = M2/M1) of the system (Peters 1989). 
Two curves (u>d and uj m in) can be drawn onto the diagram, 
where cod defines the fractional radius of the outer disk and 
Umin represents the minimum distance between the stream 
and the center of the primary star. 

We have constructed the n — q diagram (Fig. 12) for 
all the SDBs from the list of Ibanoglu et al. (2006), taking 
n and q from the catalog of Budding et al. (2004). The po- 
sitions of the gainers of S Equ and KO Aql are marked with 
big dots. The short-period systems (Port < 5-6 days) are 
located above the u>d curve; the gas stream strikes the pri- 
mary directly, resulting into an impact region and possibly 
a transient disk or annulus because the orbital separation 
of two components is too small to form a stable disk. Con- 
versely, the systems falling under the uo m in curve are long- 



period ones with a small fractional radius n and can form a 
classical accretion disk because the stream does not hit the 
primary. For the systems lying between the two curves, the 
accretion structure around the primary is very unstable and 
transient (Peters 1989, Richards & Albright 1999). S Equ 
and KO Aql belong to this group, thus they should contain 
transient, unstable disks or an accretions annulus around 
the gainers and show evidence of impact regions. 

Indeed, our Ha difference profiles suggest accretion 
structures more closely connected to the gas stream or to 
the impact region of the stream on the outer atmosphere of 
the hot star rather than to disks. The extra-absorption ob- 
served in the S Equ spectra follows the motion of the hotter 
component and is close to the star-stream impact site, as al- 
ready found by Richards (2004) who also found evidence of a 
circum-primary equatorial 'bulge' centered on the velocity of 
the primary star. The emission component (B) found by us 
displays instead a different behaviour, being likely located 
between the two stars and possibly connected to the gas 
stream itself or to the interaction between the stream and 
a transient disk. This seems to indicate a variability of the 
mass-transfer process in this system, in agreement with the 
observation of Algol systems that alternate between stream- 
like and disk-like states (Richards 2004). The chromospheric 
emission of the cool component of S Equ is also detected 
in our spectra, as already found by Richards (2001). We 
have detected similar features in the residual Ha spectra 
of KO Aql, with the exception of the excess B emission be- 
tween the two components. In addition, the extra-absorption 
displays larger velocity variations compared to S Equ. 

In Algol systems, the rotational velocity of the mass- 
gaining components is an indicator of mass transfer, which 
acts as a spin-up mechanism (Olson & Etzel 1994, Mukher- 
jee et al. 1996). Thus, when one sees a primary component 
rotating faster than expected from tidally-induced synchro- 
nism, one can argue that mass transfer has been taking 
place. However, this is not a direct proof for an ongoing mass 
flow. By the cross-correlation of the spectra of S Equ and 
KO Aql with slowly-rotating templates, we have measured, 
with a rather good accuracy (2-4 kms -1 ), the rotational ve- 
locities of the hotter components of both systems and found 



that these components are rotating about 30% faster than 
the synchronized velocity which is likely the result of the 
accretion of gas with high angular momentum. 

The mass transfer rates in Algols can be estimated as- 
suming a conservative mass and angular momentum trans- 
fer, i.e. without losses from the system. In such hypothesis, 
a monotonic increase of the orbital period is expected. An 
increasing P or b has been indeed found for S Equ (Soydugan 
et al. 2003) and KO Aql (Soydugan & Soydugan 2007). The 
analysis of the orbital period changes has allowed us to calcu- 
late the mass transfer rates dM/dt ~ 4.0 x 10~ 8 Mq j/r -1 for 

5 Equ and dM/dt ~ 2.6 x W' 7 Mqv^ 1 for KO Aql. Consid- 
ering that Algols are mostly in the slow phase of mass trans- 
fer with rates dM/dt ~ 10~ n - W~ 7 MQyr' 1 (Richards 

6 Albright 1999), the mass transfer rates for S Equ and 
KO Aql are among the highest occurring in these systems, 
thus explaining the strong Ha extra-absorption and/or emis- 
sion features shown by these systems. 

8 CONCLUSIONS 

We have observed spectroscopically the classical Algol-type 
binaries S Equ and KO Aql with the aim of determining 
the physical parameters of the components and to study the 
mass transfer and accretion structures around the hotter 
components by analyzing Ha difference profiles. 

Our detailed spectroscopic analysis has yielded the first 
precise determinations of the absolute parameters of the sys- 
tems and, in particular, the stellar masses, thanks to the 
first radial velocity curves of the cooler secondary compo- 
nents. We found that the hotter, mass-gaining components 
of both systems are underluminous compared to normal 
main-sequence stars of the same mass. We show that this 
behaviour can be explained by the Roche-lobe overflow. 

We have also found clear evidence for mass transfer, ac- 
cretion structures, and chromospheric activity of late-type 
components in our spectra. In particular, we have detected 
in S Equ a broad Ha extra-emission component that is ob- 
served red-shifted at phases around 0.25 and blue-shifted at 
phases near 0.75. This feature is not observed during the 
eclipses and its RV variation is consistent with an accretion 
structure (mass flow or interaction between the stream and 
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a transient disk) located between the two stars at about 
3x10 km from the barycenter toward the cool component. 
Similar features were also observed by Vesper et al. (2001) in 
some short-period Algol binaries such as SS Cet, U CrB, (3 
Per, and S Equ itself. However, we did not observe a similar 
behaviour in KO Aql. 

Our Ha spectra also display, for both S Equ and 
KO Aql, extra-absorption features which get stronger 
around phase and whose radial velocities are consistent 
with those of the hotter components. They are very likely 
produced by the impact region between the stream and the 
primary. Analogous features were already observed in other 
short-period Algols such as f3 Per (Richards 1993), SW Cyg 
(Hunt 1997) and RW Mon (Vesper & Honeycutt 1999). 

In the difference profiles of both systems we have also 
detected emission lines which closely follow the radial ve- 
locity curves of the late type secondaries and, as such, they 
very likely arise from their chromospheres. 
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